The performance of ensemble-based free energy protocols in computing binding affinities to ROS1 kinase

Optimization of binding affinities for compounds to their target protein is a primary objective in drug discovery. Herein we report on a collaborative study that evaluates a set of compounds binding to ROS1 kinase. We use ESMACS (enhanced sampling of molecular dynamics with approximation of continuum solvent) and TIES (thermodynamic integration with enhanced sampling) protocols to rank the binding free energies. The predicted binding free energies from ESMACS simulations show good correlations with experimental data for subsets of the compounds. Consistent binding free energy differences are generated for TIES and ESMACS. Although an unexplained overestimation exists, we obtain excellent statistical rankings across the set of compounds from the TIES protocol, with a Pearson correlation coefficient of 0.90 between calculated and experimental activities

The size matters? A computational tool to design bivalent ligands

Bivalent ligands are increasingly important such as for targeting G protein-coupled receptor (GPCR) dimers or proteolysis targeting chimeras (PROTACs). They contain two pharmacophoric units that simultaneously bind in their corresponding binding sites, connected with a spacer chain. Here, we report a molecular modelling tool that links the pharmacophore units via the shortest pathway along the receptors van der Waals surface and then scores the solutions providing prioritization for the design of new bivalent ligands. Bivalent ligands of known dimers of GPCRs, PROTACs and a model bivalent antibody/antigen system were analysed. The tool could rapidly assess the preferred linker length for the different systems and recapitulated the best reported results. In the case of GPCR dimers the results suggest that in some cases these ligands might bind to a secondary binding site at the extracellular entrance (vestibule or allosteric site) instead of the orthosteric binding site

Selective inhibitors of the PSEN1–gamma-secretase complex

Clinical development of γ-secretases, a family of intramembrane cleaving proteases, as therapeutic targets for a variety of disorders including cancer and Alzheimer’s disease was aborted because of serious mechanism-based side effects in the phase III trials of unselective inhibitors. Selective inhibition of specific γ-secretase complexes, containing either PSEN1 or PSEN2 as the catalytic subunit and APH1A or APH1B as supporting subunits, does provide a feasible therapeutic window in preclinical models of these disorders. We explore here the pharmacophoric features required for PSEN1 versus PSEN2 selective inhibition. We synthesized a series of brain penetrant 2-azabicyclo[2,2,2]octane sulfonamides and identified a compound with low nanomolar potency and high selectivity (>250-fold) toward the PSEN1–APH1B subcomplex versus PSEN2 subcomplexes. We used modeling and site-directed mutagenesis to identify critical amino acids along the entry part of this inhibitor into the catalytic site of PSEN1. Specific targeting one of the different γ-secretase complexes might provide safer drugs in the future

Best practices for constructing, preparing, and evaluating protein-ligand binding affinity benchmarks

Free energy calculations are rapidly becoming indispensable in structure-enabled drug discovery programs. As new methods, force fields, and implementations are developed, assessing their expected accuracy on real-world systems (benchmarking) becomes critical to provide users with an assessment of the accuracy expected when these methods are applied within their domain of applicability, and developers with a way to assess the expected impact of new methodologies. These assessments require construction of a benchmark - a set of well-prepared, high quality systems with corresponding experimental measurements designed to ensure the resulting calculations provide a realistic assessment of expected performance when these methods are deployed within their domains of applicability. To date, the community has not yet adopted a common standardized benchmark, and existing benchmark reports suffer from a myriad of issues, including poor data quality, limited statistical power, and statistically deficient analyses, all of which can conspire to produce benchmarks that are poorly predictive of real-world performance. Here, we address these issues by presenting guidelines for (1) curating experimental data to develop meaningful benchmark sets, (2) preparing benchmark inputs according to best practices to facilitate widespread adoption, and (3) analysis of the resulting predictions to enable statistically meaningful comparisons among methods and force fields

Selective inhibitors of the PSEN1-gamma-secretase complex

Clinical development of Y-secretases, a family of intramembrane cleaving proteases, as therapeutic targets for a variety of disorders including cancer and Alzheimer’s disease was aborted because of serious mechanism-based side effects in the phase III trials of unselective inhibitors. Selective inhibition of specific Y-secretase complexes, containing either PSEN1 or PSEN2 as the catalytic subunit and APH1A or APH1B as supporting subunits, does provide a feasible therapeutic window in preclinical models of these disorders. We explore here the pharmacophoric features required for PSEN1 versus PSEN2 selective inhibition. We synthesized a series of brain penetrant 2-azabicyclo[2,2,2]octane sulfonamides and identified a compound with low nanomolar potency and high selectivity (>250-fold) toward the PSEN1–APH1B subcomplex versus PSEN2 subcomplexes. We used modeling and site-directed mutagenesis to identify critical amino acids along the entry part of this inhibitor into the catalytic site of PSEN1. Specific targeting one of the different Y-secretase complexes might provide safer drugs in the future.The work was supported by an AIO-project (no. HBC.2016.0884). This project received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 Research and Innovation Programme (grant agreement no. ERC-834682 CELLPHASE_AD). This work was supported by the Flanders Institute for Biotechnology (VIB vzw), a Methusalem grant from KU Leuven and the Flemish Government, the Fonds voor Wetenschappelijk Onderzoek, KU Leuven, The Queen Elisabeth Medical Foundation for Neurosciences, the Opening the Future campaign of the Leuven Universitair Fonds, the Belgian Alzheimer Research Foundation (SAO-FRA), and the Alzheimer’s Association USA.Peer ReviewedPostprint (published version

Selective inhibition of intestinal guanosine 3,5-cyclic monophosphate signaling by small-molecule protein kinase inhibitors

The guanosine 3,5-cyclic monophosphate (cGMP)-dependent protein kinase II (cGKII) serine/threonine kinase relays signaling through guanylyl cyclase C (GCC) to control intestinal fluid homeostasis. Here, we report the discovery of small-molecule inhibitors of cGKII. These inhibitors were imidazole-aminopyrimidines, which blocked recombinant human cGKII at submicromolar concentrations but exhibited comparatively little activity toward the phylogenetically related protein kinases cGKI and cAMP-dependent protein kinase (PKA). Whereas aminopyrimidyl motifs are common in protein kinase inhibitors, molecular modeling of these imidazole-aminopyrimidines in the ATP-binding pocket of cGKII indicated an unconventional binding mode that directs their amine substituent into a narrow pocket delineated by hydrophobic residues of the hinge and the C-helix. Crucially, this set of residues included the Leu-530 gatekeeper, which is not conserved in cGKI and PKA. In intestinal organoids, these compounds blocked cGKII-dependent phosphorylation of the vasodilator-stimulated phosphoprotein (VASP). In mouse small intestinal tissue, cGKII inhibition significantly attenuated the anion secretory response provoked by the GCC-activating bacterial heat-stable toxin (STa), a frequent cause of infectious secretory diarrhea. In contrast, both PKA-dependent VASP phosphorylation and intestinal anion secretion were unaffected by treatment with these compounds, whereas experiments with T84 cells indicated that they weakly inhibit the activity of cAMP-hydrolyzing phosphodiesterases. As these protein kinase inhibitors are the first to display selective inhibition of cGKII, they may expedite research on cGMP signaling and may aid future development of therapeutics for managing diarrheal disease and other pathogenic syndromes that involve cGKII